INCONTRI SULLA FISICA DELLE ALTE ENERGIE 



Torino, April 14-16, 2004 



Status and prospects of charm physics: a few 

considerations 

S. Malvezzi^) 
1. INFN, Milano, Italy 

Abstract — The goals of quark flavor physics are to test the consistency of the Standard Model (SM) 
description of quark mixing and CP violation, to search for evidence of New Physics, and to select 
between New Physics scenarios that might be initially uncovered at the LHC. This will require a range 
of measurements at the percent level or better in the flavor-changing sector of the SM. Charm repre- 
sents a unique candidate for these studies. Progress in charm physics has been prodigious over the 
last twenty years; it comes both from fixed-target and collider (e+e^ and pp) experiments, and can 
guide us toward future investigations. To fully exploit the potentiality of flavor physics in indicating 
New Physics, non-perturbative strong-interaction dynamics has to be dealt with. Complications in the 
determination of CP phases due to strong final-state interaction (FSl) have been uncovered; studies 
of charm-meson decays, mainly through Dalitz-plot analyses, have experimentally confirmed the rel- 
evance of FSI phases. Analogously, determination of Cabibbo-suppressed CKM matrix elements and 
DD mixing parameters will require an understanding of strong-interaction effects among the light 
hadrons produced in heavy-meson decays. The same complications will affect the beauty sector and 
will need to be kept under control to search for New Physics effects. The role of charm in the search 
for New Physics, the lessons learnt so far and warnings for the future will be discussed. 

1 Beyond the Standard Model; the clue from charm 

The success of the Standard Model in describing the experimental information to date suggests that possible devi- 
ations have to be pursued either at high-energy scales or as small effects in low-energy variables. In the Standard 
Model, CP violation (CPV) and flavor-changing neutral currents (FCNC) are expected to be small in charm decay, 
while lepton-family number violation and lepton-number violation are forbidden. Anomalously large rates for 
rare decays would imply non-Standard Model tree-level FCNC diagrams or non-Standard Model contributions to 
higher-order loop diagrams, which make them sensitive to high-mass gauge bosons and fermions and allow for 
probing particle states and mass scales not directly accessible [1]. The study of FCNC has been mainly dedicated 
to transitions such as s ^ dl^l^ , s dvv, b ^ s-f and b sl^l^, as well as phenomena such as i^'^-A''^ and 
50_gO ijyxing. The analogous FCNC processes in the charm sector have been investigated less. The Standard- 
Model expectations for both D^-D^ mixing and FCNC are very small, while extensions of the Standard Model 
may enhance the contributions, which can then be orders of magnitude larger. In the Standard Model, the D sys- 
tem is not as sensitive to CP as the K and B mesons. In this case too, the small predicted effects could leave 
open a window to the observation of New-Physics effects. Charm decays may be the only window onto such New 
Physics since it is possible that the mechanism responsible will only couple to up-type quarks. More specifically, 
non-Standard-Model forces might exhibit very different patterns for the up and down classes of quarks [2]. Charm 
decays are the only up-type decays that afford a probe of such physics: non-strange light-flavor hadrons do not 
allow for oscillations and top-flavored hadrons do not even form in a practical way. In Table 1 experiments that 
have already published results on charm mixing, rare decays and CPV searches are listed, along with their charac- 
teristics and a reference sample. In Table 2 the expected data sets for existing experiments and proposed facihties 
are reported [3]. 



Table 1: A comparison of the LEP experiments, CDF, E791, FOCUS, CLEO, BaBar and Belle. K-'k+ is the 
number of reconstructed K~-k^ used in the pubUshed measurements, at denotes the proper time resolution 
for charm hadrons. 
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Table 2: Expected K tt'^ samples for the existing experiments and proposed facilities. K 7r+ is the 

number of of reconstructed D''^ K^ir^ in the full data set. 



Expcrimcnl 


Full data 
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K 7r+ rcconslructcd 


BaBar 
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-1 


6.6 X 10^' 


Belle 


500 fb- 


-1 


6.6 X 10'^ 


CDF 


4.4 fb" 


■ 1 


30 X 10^ 


CLEO-c 


3.0fb- 
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5.5 X IQS 


BES III 


30.0 fb- 


-1 


5.5 X 10^ 


BTeV 






~ 6 X 10^/10^ 


SuperBaBar 


10.0 ab- 


-1 


1.3 X 10^/10'' s 


SuperKEKB 


2.0 ab- 


-1 


2.5 X 10710^ s 



1.1 Mixing 

£)0 _ f)0 juixing in the Standard Model occurs because the two eigenstates and D° are not mass eigenstates. 
The probabiUty that a meson produced at t = decays as a at time t is given by 

P(D0 ^ D°) = ie-i^i* {1 - 2e^ cos Ami + e"^"} (1) 

where Am = m2 — mi, AF = — Fi are the mass and width differences of the mass eigenstates and F = ^'^2^^ . 
In the limit Am, AF << F, mixing effects are parametrized by two dimensionless amplitudes x = and 
y = The experimental measurement is performed by analyzing wrong sign semileptonic and hadronic decays; 
the presence of doubly Cabibbo-suppressed decays complicates the measurement in the hadronic sector. Assuming 
CP conservation, the wrong sign to right sign decay rate is 



roit) = e-r* (rdcs + VRdcs y'rt + ^{x'^ + y'^)rH^ 



(2) 



where y' = y cos S — x sin S and x' = x cos 6 + y sin S, S being the strong phase difference between the Cabibbo- 
favored and Doubly Cabibbo-suppressed Decay (DCS). The three terms come from DCS decays, interference, and 
mixing. In the semileptonic mixing only the last term matters. The up-to-date experimental situation is summarized 
in Fig. 1 [4, 5]. In the semileptonic sector, FOCUS has measured, preliminarily, r/j in the mode D'^ K'^ ii~v^ 
to be less than 0.131% @ 95%C.L.; CLEO-c expects td less than @ 95%C.L.; SuperBaBar estimates 

that, with 10 ab~^, the sensitivity would be < 5 x 10~^. With 10^ reconstructed charm decays BTeV expects to 
reach a sensitivity of 1-2 x IQ- ''. An observation of mixing at the level of 10^'' — IQ-"'', which could be probed by 
a high-sensitivity experiment, would be a signal of New Physics. In Fig. 1 measurements of the hfetime difference 
between decays to CP-even and CP-odd final states are also shown. If CP violation in neutral Z)-meson decays 
is negligible we can write 

_ r(CP even) - r(CP odd) r(£)° ^ K+R-) _ 

~ F(CP even) + r(CP odd) ~ T{D^ ^ K-tt+) ^ ' 

The theoretical prediction spread of the mixing parameters is wild, as shown in the right side of Fig. 1; indepen- 
dently of New Physics discovery, accurate experimental measurements are necessary. 



1.2 CP Violation 

CP violation can occur in charm decays via the interference of two weak decay amplitudes to the same final 
state. Indirect CP violation is mediated by D^-D^ mixing. The interference is the largest if the two amplitudes 
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Figure 1 : Mixing scenario. 

are roughly equal. However, since mixing is expected to be small, indirect CP violation is not expected to be a 
big effect. In absence of mixing, a decay mode (e.g., a Cabibbo-suppressed decay) that has two weak amplitudes 
contributing to the same final state can exhibit direct CP violation. Final State Interactions (FSI) can induce a phase 
shift between the two weak amplitudes, leading to a decay-rate asymmetry between a charm meson decay and its 
CP conjugate: r(D — > /) 7^ r{D /). FSI are substantial in charm decay, and the two weak amplitudes can be 
similar in size. Asymmetries as large as 0.1-1% are possible in the Standard Model. In Table 3 the experimental 
scenario is reported. 



Table 3: Measurements of CP asymmetries. 
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-0.017 ±0.042 












FOCUS 
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D°- 
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-0 


001 ± 0.022 ±0.015 





020 ±0.012 ±0.006 


D°- 


-> 7r^7r+ 





048 ±0.039 ±0.025 





030 ±0.013 ±0.006 


D+ - 


K-K+TT+ 





006 ±0.011 ±0.005 






D+ - 




-0 


016 ±0.015 ±0.009 






D+ - 







071 ±0.061 ±0.012 







Expectations come from CLEO-c, which aims at exploiting the quantum coherence of the D'^-D'^ pair produced. 
The process e+e^ ijj" D^D'^ produces a CP+ eigenstate in the first step, since ij}" has J^'~^ ~ 1 . 
Consider the case where both the and the decay into CP eigenstates. Then the decays i/i f^fX or 
/!_/!_ are forbidden, where f+/ f- denotes a CP+ / CP- eigenstates. This is because CP{fj.fj^) = (—1)' = —1 
forthel=l ip . Thus, observation of a final state such as (/'ir+A'^)(7r+7r^) constitutes evidenceofCPviolation. Of 
course this method requires a complete understanding of the Initial State Radiation (ISR) effects. CP asymmetries 
measured so far in the charm sector are consistent with zero within the errors. The high statistics soon available 
from experiments such as Belle, BaBar and CLEO-c would be able to probe CP asymmetries at the 10^'^ level. 
BTeV and will extend the sensitivity below the lO^'^ level. 
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1.2.1 CP investigation via Dalitz-plot analysis 



Dalitz-plot analysis has recently emerged as a unique tool to investigate CPV effects. Amplitude analysis gives 
the full observation of a three-body decay, providing access to the coefficients and phases of the various channels. 
CP effects are intimately connected to phase phenomena. The comparison between the relative phases of the 
CP-conjugate states for the various resonant channels can address CPV effects. FOCUS analyzed the _D+ 
K~K^Tr^ channel [6] and CLEO the Kstt^tt^ channel [7]. In general each measured phase in the 

Dalitz plot can be interpreted as a sum of two components, one CP conserving, the other CP violating. Under 



CP conjugation the first does not switch sign, while the second does. The 



K K^TT^ channel is a good 



candidate to observe CP effects. This decay is Cabibbo suppressed and thus contains two diagrams contributing to 
the final state, spectators and penguins; it is also dominated by strong effects, as discussed above, and thus satisfies 
the two necessary requirements for detecting a CP asymmetry. FOCUS [6] has performed a new Dalitz-plot study, 
which proceeds through a complete Dalitz analysis of the / D~ split samples. The preliminary results on 
/ coefficients and phases, shown in Fig. 2, do not point to any CPV manifestation. The relative phases 
are measured with an accuracy of about 5-10°. The errors are affected by a large correlation among the various 
resonant components necessary to fit the Dalitz plot. A more robust formalism, such as the K-matrix approach, 
which will be discussed later in this paper, might allow for a good fit with fewer resonances, thus ameliorating the 
accuracy of the phase measurements. At any rate a factor of 100 more statistics, as hopefully available in BTeV, 
will allow for reaching sensitivities of a degree. 
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Figure 2: Amplitude coefficient and relative phases (in degrees) of all the resonant contributions of the 
K^TT^ decay: the circles are for the full sample, upward triangles for and downward triangles for 



1.3 Rare decays 
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Figure 3: Measurements of rare decays. 

The search for rare and forbidden decays is enticing since Standard Model predictions tend to be beyond the reach 
of current experiments, and a signal is a sign of unexpected physics. Standard Model predictions of rare decays 
are dominated by long-range effects, which are notoriously difficult to calculate and experimental measurements 
are thus crucial. The current experimental scenario concerning rare decay branching ratios is reported in Fig. 3. 
FOCUS [8] has improved many of the previous limits of a factor from 2 to 14. The strength of the FOCUS fixed 
target experiment is the excellent vertexing, which was used to require that the two candidate leptons form a good, 
well separated vertex from the primary. CDF and DO have efficient di-lepton triggers and future measurements are 
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promising. The reported CDF result is obtained only with a statistics of 65pb ^. A new limit of 2 x 10 ^ has 
been recently provided by the Hera-B experiment [9]. CLEO-c sensitivity is estimated at the level of 10^^. 



2 Standard Model measurements and lessons for the future 

The variety of measurements above briefly discussed represents an excellent opportunity for searching New Physics 
in the charm sector. Effects beyond the Standard Model are, as already pointed out, expected to be small, making 
their investigations really challenging. Possible tiny inconsistencies and deviations from expected behavior can 
be interpreted as sign of New Physics only if the Standard Model phenomenology is correctly parametrized in 
the analysis procedure. The main lesson coming from the mature field of charm is that the advantages of high 
statistics to interpret the Heavy Flavor decay dynamics will vanish in the absence of a strategy to control strong 
effects among particles involved in weak decay processes. A few examples from the FOCUS experiment will be 
discussed. 



2.1 Unexpected new results in the semilep tonic sector 

Semileptonic decays have always been considered the best candidates to study the charm phenomenology. Decay 
rates can be calculated from first principles, e.g, Feynman diagrams; due to the presence of only one hadron in 
the final state we do not have to worry much about final state interaction complications; the hadronic part of the 
decay can be contained in proper form factors, which, on their turn can be predicted by theory such as HQET, 
LGT and Quark models. However recent FOCUS studies have shown that decays in the semileptonic sector 
also reveal the presence of quantum mechanical effects and hadronic complications have to be dealt with in the 
analysis. Of particular relevance is the study of I?+ X^vr+zi+i/ [10]. The FOCUS mass spectrum appears, as 
expected, dominated by A'*°(892); yet an unexpected forward-backward asymmetry was exposed in the cosOy 
variable, defined in Fig. 4. The asymmetry is striking below the pole and essentially absent above the pole: it thus 
suggests quantum-mechanical interference between the Breit-Wigner amplitude describing the K*'^{892) and a 
broad or nearly constant S'-wave amplitude. FOCUS has tried to describe this behavior of the decay distribution 
for K^TT^fi'^u through a simple model. It has been written in terms of the three helicity-basis form 

factors: Ha, and the S'-wave amplitude modeled as a constant, with modulus A and phase S has been 
added. Angular momentum conservation restricts the S'-wave contribution to the Hq piece that describes the 
amplitude for having the virtual in a zero helicity state relative to its momentum vector in the D+ rest-frame. 




Figure 4: The definition of the kinematic variables for the £)+ — > K*^ fi^v and cos Oy distribution. 

The shape of the cos Oy term versus wlk-k turns out to be a strong function of the interfering S-wave amplitude 
phase 5. Figure 5 shows how the cos By distribution (properly weighted) is consistent with a constant S-wave 
amplitude of the form 0.36 exp(j7r/4) (GeV)^^, although, at this stage of the analysis, the solution is not unique; 
alternative modelings of the S-wave amplitude, as shown in the same Fig. 5, could also fit the data. The measured 
phase of 7r/4 is consistent with that found by LASS from the Kit phase-shift analysis. The hypothesis of the 
broad scalar resonance k is disfavored since it would imply a phase-shift of approximately 90"^ with respect to 
A'*°(892), unlikely to manifest in the semileptonic sector where FSI are not expected to be wild [11]. A more 
precise determination can be obtained through a complete analysis of the form factor ratios and r2 [12]. FOCUS 
has found r„ = 1.504 ± 0.057 ± 0.039 and ra = 0.875 ± 0.049 ± 0.064. For the S-wave we obtained an 
amplitude modulus of A 0.330 ± 0.022 ±0.015 GeV"! and a phase (5 = 0.68 ± 0.07 ± 0.05 rad, in reasonable 
agreement with the very informal previous estimate of A — 0.36 exp(i7r/4) GeV^^. The inclusion of the S-wave 
amplitude has dramatically improved the quality of the form factor fits. A measurement of the branching ratio 
Y{^l^+ ^K- Ji t!+-) = 0.602 ± 0.01 (stat.) ± 0.021 (syst.) has then been performed including interference with this 
broad scalar resonance. To complete the semileptonic study, FOCUS has performed an analysis of the Df 
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Figure 5: Asymmetry distribution in K-k invariant mass. The dashed Une is the MC with no S'-wave ampHtude, the 
soHd line is the MC with an S'-wave ampHtude of 0.36 exp(j7r/4) (GeV)^^. Different modehngs of the S'-wave 
are also shown. 



<^li^v final state; the ^^ft 



r(r>j^</.7r+) 



has been measured to be 0.54 ± 0.033(stat.) ± 0.048(syst.) [13], and a new 



measurement of form factor ratios has been recently published [14]. 



2.2 Unexpected complications in hadronic decay dynamics: and Z)+ ^ tt+tt 7r+ 

Charm-meson decay dynamics has been extensively studied in the last decade. The analysis of the three-body 
final state by fitting Dalitz plots has proved to be a powerful tool for investigating effects of resonant substruc- 
ture, interference patterns, and final state interactions in the charm sector. An example of Dalitz-plot analysis to 
study CP-violation effects has been already discussed in a section of this paper. The isobar formalism, which has 
traditionally been applied to charm amplitude analyses, represents the decay amplitude as a sum of relativistic 
Breit-Wigner propagators multiplied by form factors plus a term describing the angular distribution of the two 
body decay of each intermediate state of a given spin. Many amplitude analyses require detailed knowledge of 
the light-meson sector. In particular, the need to model intermediate scalar particles contributing to the charm 
meson decays into three-body hadronic channels has caused experimentalists of the field to question the validity 
of the Breit-Wigner approximation for the description of the relevant scalar resonances [15, 16, 17]. Resonances 
are associated with poles of the S-matrix in the complex energy plane. The position of the pole in the complex 
energy plane provides the fundamental, model-independent, process-independent resonance description. A simple 
Breit-Wigner amplitude corresponds to the most elementary type of extrapolation from the physical region to an 
unphysical-sheet pole. In the case of a narrow, isolated resonance, there is a close connection between the position 
of the pole on the unphysical sheet and the peak we observe in experiments at real values of the energy. However, 
when a resonance is broad and overlaps with other resonances, then this connection is lost. The Breit-Wigner 
parameters measured on the real axis (mass and width) can be connected to the pole-positions in the complex en- 
ergy plane only through models of analytic continuation. A formalism for studying overlapping and many channel 
resonances has been proposed long ago and is based on the K-matrix [18, 19] parametrization. This formalism, 
originating in the context of two-body scattering, can be generalized to cover the case of production of resonances 
in more complex reactions [20], with the assumption that the two-body system in the final state is an isolated one 
and that the two particles do not simultaneously interact with the rest of the final state in the production process 
[19]. The K-matrix approach allows for including the positions of the poles in the complex plane directly in the 
analysis, incorporating in the charm analysis the results from light spectroscopy experiments [21, 22]. In addition, 
the K-matrix formalism provides a direct way of imposing the two-body unitarity constraint which is not explicitly 
guaranteed in the simple isobar model. Minor unitarity violations are expected for narrow, isolated resonances but 
more severe ones exist for broad, overlapping states. The validity of the assumed quasi two-body nature of the 
process of the K-matrix approach can only be verified by a direct comparison of the model predictions with data. 
In particular, the failure to reproduce three-body-decay features would be a strong indication of the presence of 
neglected three-body effects. 



2.3 The isobar formalism 

The formalism traditionally applied to three-body charm decays relies on the so-called isobar model. A resonant 
amplitude for a quasi-two-body channel, of the type 

D ^ r + c 
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is interpreted a la Feynman. For the decay D 

TT, 



TTTTTT of Fig. 6, a D 




TT current with form factor Fjj interacts 

+ 



(P7r+ -P7r-)i^ 



Figure 6: The tttttt decay diagram. 

with a di-pion current with form factor Fr through an unstable propagator with an imaginary width contribution in 
the propagator mass. Each resonant decay function is thus, 



A = FoFr X \c\'\a\'Pj {cos ei J x BW{mab) 



(5) 



i.e., the product of two vertex form factors (Blatt-Weisskopf momentum-dependent factors), a Legendre poly- 
nomial of order J representing the angular decay wave function, and a relativistic Breit-Wigner (BW). In this 
approach, already applied in the previous analyses of the same channels [23], the total amplitude (Eq. 6) is as- 
sumed to consist of a constant term describing the direct non-resonant three-body decay and a sum of functions 
(Eq. 5) representing intermediate two-body resonances. 



A{D) = aoe 



(6) 



2.4 The K-matrix formalism 

For a well-defined wave of specific isospin and spin //, characterized by narrow and isolated resonances the propa- 
gator is, as anticipated, of the simple BW form. In contrast, when the specific wave // is characterized by large and 
heavily overlapping resonances, just as the scalars, the propagation is no longer dominated by a single resonance, 
but is the result of complicated interplay among the various resonances. In this case, it can be demonstrated on 
very general grounds that the propagator may be written in the context of the K-matrix approach as 



{I-iK-pY 



(7) 



where K is the matrix for the scattering of particle a and b ( Eq. 4) and p is the phase-space matrix. In this picture, 
the production process is viewed as consisting of an initial preparation of several states, which then propagate via 
the term (/ — iKp)~^ into the final state. In particular, the three-pion final state can be fed by an initial formation 
of {Tnr)n, {KK)w, {'qrf)'K, {'q'q')-K and multi-meson states (mainly four-pion states at ^/s < 1.6 GeV). While the 
need for a K-matrix parametrization, or in general for a more accurate description than the isobar model, might be 
questionable for the vector and tensor amplitudes, since the resonances are relatively narrow and well isolated, this 
parametrization is needed for the correct treatment of scalar amphtudes. Indeed the tttt scalar resonances are large 
and overlap each other in such a way that it is impossible to single out the effect of any one of them on the real axis. 
In order to write down the propagator, we need the scattering matrix. To perform a meaningful fit to D mesons to 
three-pion data, a full description of the scalar resonances in the relevant energy range, updated to the most recent 
measurements in this sector is needed. At the present time the only self-consistent description of S'-wave isoscalar 
scattering is that given in the K-matrix representation by Anisovich and Sarantsev in [22] through a global fit of the 
available scattering data from the tttt threshold up to 1900 MeV. FOCUS has performed the first fit to charm data 
with the K-matrix formalism [24] in the Z)+ and 2?+ — > 7r+7r~7r+ channels. The K-matrix used is that of [22]: 



(a) (a) 



xscatt 
lij 



1 GeV^ - 8%'=''^* 



S - Sn 



s - SAml/2 
{s - SAo)(l - Sao)" 



(8) 



The factor f;!"^ is the coupling constant of the K-matrix pole a to meson channel i; the parameters 
gscatt (jgsQj-ibg a smooth part of the K-matrix elements; the factor (^^g^^J^^j^l^.^^) suppresses a false kinematical 
singularity in the physical region near the tttt threshold (Adler zero). The K-matrix values of [22] generate a phys- 
ical T-matrix, T = {I — ip ■ K)~^K, which describes the scattering in the (00)"'""'"-wave with five poles, whose 



scatt j^jjj 
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masses, half-widths, in GeV are (1.019,0.038), (1.306,0.167), (1.470,0.960), (1.489,0.058) and (1.749,0.165). 
The K-matrix formalism, originated in the context of the two-body scattering, can be generalized to deal with 
formation of resonances in more complex reactions, through the P-vector [20] approach. The decay amplitude for 
the D meson into three-pion final state, where tt+tt^ are in a {IJ^'~^ = 00++)-wave, can thus written as 



AiD {7T+7r-)oo^+7T+) ^ 



E 



Pa9 



(a) 



/•prod 



IGeV^ 



prod 



prod 



s ~ SAml/2 
(s - SAo)(l - SAo) ' 



(9) 



where f3a is the coupling to the pole a in the 'initial' production process, /f ™^ and Sq"^"" are the P-vector slowly 



prod 



varying (SVP) parameters. In the end, the complete decay amplitude of the D meson into three-pion final state is 
[24]: 

A{D) = aoe'^° + ^ ate'^'Ai + Fi (10) 



where the index i now runs only over the vector and tensor resonances, which can be safely treated as simple 
Breit-Wigner's (see Eq. 5). In the fit to the data, the K-matrix parameters are fixed to the values of [22], which 
consistently reproduce measured S'-wave isoscalar scattering. The free parameters are those peculiar to the P- 



vector, i.e., /3q, ff^°'^ and Sq'°°, and those in the remaining isobar part of the amplitude, and Si. 



prod 
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Figure 7: Signal and side-band regions of the three-pion invariant-mass distribution for a) Df and b) Dalitz- 
plot analysis respectively. 
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Figure 8: a) D+ and b) D+ DaHtz plots. 
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The three-pion selected samples (Fig. 7) consist of 1527± 51 and 1475 ± 50 events for the and respectively. 
The Dalitz-plot (Fig. 8) analyses are performed on yields within ±2(t of the fitted mass value. 

2.4.1 Results for the £)+ — > 7r+7r~7r+ decay 

The general procedure adopted for the fits consists of several successive steps in order to eliminate contributions 
whose effects on fits are marginal. Initially all the well established, non-scalar resonances decaying to tt+tt^ with a 
sizeable branching ratio are considered. Contributions are removed if their amplitude coefficients, ai of Eq. 10, are 
less than 2 a significant and the fit confidence level increases due to the decreased number of degrees of freedom 
in the fit. The P-vector initial form includes the complete set of K-matrix poles and slowly varying function (SVP) 
as given in reference [22]; j3a as well as the f^™'^ terms of Eq. 9 are removed with the same criteria. The fit 
confidence levels (C.L.) are evaluated with a estimator over a Dalitz plot with bin size adaptively chosen to 
maintain a minimum number of events in each bin. Once the minimal set of parameters is reached, addition of 
each single contribution previously eliminated is reinstated to verify that the C.L. does not improve. The resulting 
fit fractions ^\ phases and amplitude coefficients are quoted in Table 4. Both the three-body non-resonant and 
p°(770)7r+ components are not required by the fit. This result is to be compared with that obtained with the simple 
isobar model, which requires a non-resonant component of about 25% to get a decent fit to the data [25]. This 
component, which crosses the Dalitz plot uniformly, seems to compensate, with its interference with the other 
contributions, for the inability of the model to properly describe some non-trivial resonant features not properly 
accounted for in the model. In this way the potentiality of the Dalitz-plot analysis to gauge the level of the 
annihilation contribution in the charm hadronic decays is limited. An additional difficulty with the isobar model is 
the general poor knowledge of scalar resonances: the measurements reported in the PDG are dispersed over a wide 
range of values and can not be used as input parameters of charm decay amplitudes. Masses and widths of the 
corresponding Breit-Wigner forms have to be let free in the fit: the isobar model can thus be viewed as an effective 
model able to reproduce the data with a sum of effective resonances but its phenomenological interpretation has to 
be considered with caution. The entire S'-wave contribution obtained with the K-matrix formalism is represented 
by a single fit fraction since, as previously discussed, one cannot distinguish the different resonance or SVP S'- 
wave contributions on the real axis. The 13+ Dalitz projections of FOCUS data are shown in Fig. 9 superimposed 
with final fit projections. The fit C.L. is 3%. 



Table 4: Fit results from the K-matrix model for 



decay channel 


fit fraction (%) 


phase (deg) 


amplitude coefficient 


(S'-wave) 7r+ 
/2(1270) 7r+ 
/9°(1450)7r+ 


87.04 ±5.60 ±4.17 
9.74 ± 4.49 ± 2.63 
6.56 ±3.43 ±3.31 


(fixed) 
168.0±18.7±2.5 
234.9 ± 19.5 ± 13.3 


1 (fixed) 
0.165 ±0.033 ±0.032 
0.136 ±0.030 ±0.035 


Fit C.L 


3.0% 








Low ^+ ^- mass projection (GeV^ Ic^) High ^- mass projection {GeV^ Ic^) 



Figure 9: Dalitz-plot projections with the K-matrix fit superimposed. The background shape under the signal 
is also shown. 

^'The quoted fit fractions are defined as tlie ratio between the intensity for a single amplitude integrated over the Dalitz plot and that of the 
total amplitude with all the modes and interferences present. 
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2.4.2 Results for the D+ — > 7r+7r~7r+ decay 

The _D+ 7r+7r^7r+ Dalitz plot shows an excess of events at low tt+tt^ mass, which cannot be explained 
in the context of the simple isobar model with the usual mixture of well established resonances along with a 
constant, non-resonant amplitude. A new scalar resonance, the cr(600), has been previously proposed [26] to 
describe this excess. However we know that complex structure can be generated by the interplay among the 5*- 
wave resonances and the underlying non-resonant 5'-wave component that cannot be properly described in the 
context of a simple isobar model. It is therefore interesting to study this channel with the present formalism, which 
embeds all the experimental knowledge about the S'-wave tt+tt^ scattering dynamics. With the same procedure 
based on statistical significance and fit confidence level used in the Df analysis, the final set of contributions 
is reached. Beside the S'-wave component, the decay appears to be dominated by the p'^(770) plus a /2(1270) 
component. The p°(1450) was always found to have less than 2 a significance and was therefore dropped from 
the final fit. In analogy with the Df, the direct three-body non-resonant component was not necessary since the 
S VP of the S'-wave could reproduce the entire non-resonant portion of the Dalitz plot. The complete fit results are 
reported in Table 5. The Dalitz projections are shown in Fig. 10. The fit C.L. is 7.7%. 



Table 5: Fit results from the K-matrix model fit for D 



decay channel 


fit fraction (%) 


phase (deg) 


amplitude coefficient 


(S'-wave) 7r+ 
/2(127G)7r+ 
pO(770) 7r+ 


56.00 ±3.24 ±2.08 
11.74 ± 1.90 ±0.23 
30.82 ±3.14 ±2.29 


(fixed) 
-47.5 ± 18.7 ± 11.7 
-139.4 ±16.5 ±9.9 


1 (fixed) 
1.147 ±0.291 ±0.047 
1.858 ± 0.505 ±0.033 


Fit C.L. 


7.7% 








Low ^+ ^- mass projection {GeV^ Ic^) High ^+ ^- mass projection (GeV^ Ic") 

Figure 10; Dalitz -plot projections with the final fit superimposed. The background shape under the signal is 
also shown. 

The most interesting feature of these results is the fact that the better treatment of the 5-wave contribution provided 
by the K-matrix model can reproduce the low-mass tt+tt^ structure of the Dalitz plot. This suggests that any 
cr-like object in the D decay should be consistent with the same cr-like object measured in the tt+tt^ scattering. 
Additional studies with higher statistics will be required to completely understand the a puzzle. It is interesting 
to recall the close analogy between the D tttttt channel and the B ^ pn one, which is a good candidate to 
measure the angle a of the Standard Model Unitarity Triangle; the analysis of B ^ pir will proceed through 
a time-dependent Dalitz-plot analysis of the three-pion final state and will likely present similar parametrization 
compUcations. 

2.5 Interpretation of the and D+ — ^ 7r+7r~7r+ results 

The K-matrix formalism has been applied for the first time to the charm sector in the FOCUS Dalitz-plot analyses 
of the Df and 7r+7r~7r+ final states. The results are extremely encouraging since the same K-matrix 

description gives a coherent picture of both two-body scattering measurements in light-quark experiments as well 
as charm meson decay. This result was not obvious beforehand. Furthermore, the same model is able to reproduce 
features of the _D+ 7r+7r^7r+ Dalitz plot that otherwise would require an ad hoc a resonance. In addition, the 
non-resonant component of each decay seems to be described by known two-body S'-wave dynamics without the 
need to include constant amplitude contributions. The K-matrix treatment of the 5-wave component of the decay 
amplitude allows for a direct interpretation of the decay mechanism in terms of the five virtual channels considered; 
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TTTT, KK, r]r], rjr]' and 4 n. By inserting KK ^ in the decay amplitude, F, 

F={I- iKp)-'^P = (/ - iKp)-'^KK-^P = TK-^P = TQ (1 1) 

we can view the decay as consisting of an initial production of the five virtual states which then scatter via the 
physical T into the final state. The Q-vector contains the production amplitude of each virtual channel in the 
decay. The resulting picture, for both Df and £)+ decay, is that the S'-wave decay is dominated by an initial 
production of rjrj, rjrj' and KK states. Dipion production is always much smaller. This suggests that in both 
cases the S'-wave decay amplitude primarily arises from a ss contribution such as that produced by the Cabibbo- 
favored weak diagram for the Df and one of the two possible singly Cabibbo-suppressed diagrams for the D+. 
For the D~^, the ss contribution competes with a dd contribution. That the /o(980) appears as a peak in the nn 
mass distribution in £)+ decay, as it does in Dg decay, shows that for the 5-wave component the ss contribution 
dominates [16]. Comparing the relative S'-wave fit fractions that we observe for D+ and D+ reinforces this picture. 
The S-wave decay fraction for the -D+ (87%) is larger than that for the (56%). Rather than coupling to an 
5-wave dipion, the rfJ piece prefers to couple to a vector state like /9°(770) that alone accounts for ~ 30 % of 
decay. This interpretation also bears on the role of the annihilation diagram in the Df TT~^TT~Tr~^ decay. 
This study suggests that the S-wave annihilation contribution is negligible over much of the dipion mass spectrum. 
It might be interesting to search for annihilation contributions in higher spin channels, such as p°(1450)7r and 
/2(1270)7r. 

3 Conclusions 

Dedicated studies of charm physics are foreseen in the near future at existing collider facilities such as CLEO-c, 
the iJ-factories, and the Tevatron. Other facilities might also play a role, such as GSI and Super 5-factories. In 
parallel, the mature stage of the ongoing analysis is providing warnings and lessons on how to perform proper 
studies to interpret and understand Heavy Flavor phenomenology and, hopefully, to find signs of New Physics. 
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